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1. BASIC Programming 



CLINI CA L PROGRAMS WITH GAxMMA-11 

A) BASIC Programming 

BASIC is loaded under the control of the RT-11 monitor by typing 

CTRL/C 
.R BASIC <^CR^ 

BASIC is then loaded and prints: 

BASIC Vi2flB-j2f2 
* 

and awaits the specification of the optional functions available 

to the system. The user should type a <('CR/' and BASIC will 

print: 

USER FNS LOADED 
READY 

and waits for a user command. (NOTE: refer to the BASIC/RT-11 

Reference Manual, section 1.1 for alternate responses to 

loading the user functions.) 

To load and execute a GAMMA-11 clinical program, the user must 
issue the OLD command (refer to the BASIC/RT-11 Reference 
Manual, section 7.2) to erase the storage area and input the 
desired program. The user should type: 

OLD <^CR^ 
BASIC responds with: 

OLD FILE NAME _ _ 
and waits for the program name and a carriage return. 



-2- 

BASIC will respond: 

READY 
after the program is loaded. To execute the program the 
user types: 

RUN <CR^ 
and the program begins. It is recommended that prior to re- 
executing any GAMMA-11 clinical program the user re-loads the 
program with the OLD command before issuing the RUN command. 

Users who wish to edit the GAMMA-11 clinical programs or create 
their own programs are advised to read Chapter 7 of the 
BASIC/RT-11 Language Reference Manual for a complete description 
of the BASIC editting commands . 

B) Program Descriptions 

The following clinical programs are listed by their program 
names. These names should be used in conjunction with the OLD 
command in BASIC/RT-11 to load the program into memory, 
(See EXAMPLE 1) 

EXAMPLE It 
♦R BASIC 

BASIC yOlB-02 
* 

USER FNS LOADED 

READY 

OLD 

OLD FILE NAME~-CURyE 

READY 



2. EJECTS ; Ejection fraction calculation (static method) 



EJECTS - ejection fraction calculation (static method) 
The ejection fraction can be calculated by saving two matrices 
in save areas, one matrix an image of end-diastole and the 
other an end-systolic image. 

In both save areas, a region of interest (ROI) corresponding 
to the left ventricle must be included. 

The program requests the save area number and the ROI 
corresponding to end-diastole. Then the program requests the 
end-systole save area number and ROI. The program will output 
the patient index line and the ejection fraction computed. 

The method of determining the ROI of the left ventricle in 
end-diastole and the ROI in end-systole is completely up to 
the user. One such method is the technique of Kurtz, Ahnberg, 
Freed, LaFarge and Treves (to be published) : 

1) There are two methods for acquiring an end-diastole image, 
and then an end-systole image. 

a) gating the camera output to produce one end-diastole 
image and one end-systole image.* 

b) acquiring a fast (around 10 frames/sec) dynamic study. 
Save a curve over the left ventricle and call the BASIC 
program MARK. Pick the end-diastole frames (peaks) and 

*suGh as a Brattle gate or a Riverside Bio-Engineering gate. 



the end-systole frames (valleys) from the curve displayed. 
The frame numibers are printed on the output device. Add 
the 1st column frame numbers to create an end-diastole 
image. The end-systole image is the summation of the 2^^ 
column. (Figure 6) 

2) Display the end-diastole image, lower the upper threshold. 

(recommended setting is between 45 and 65%) (e.g. UT55) 

3) Smooth the data (S9) 

4) Enter irregular region of interest (IR) , outline the left 
ventricle and save the matrix (SMr,) . (Figure 1) 

5) Display the end-systole image, lower the upper threshold 
(recommended setting is around 5% lower than threshold in 
step 2) (e.g. UT48) and repeat steps 3 and 4. (Figure 2) 

6) Call program EJECTS. 




FIGURE 1. ROI defined over end-diastole left 
ventricle 




FIGURE 2. ROI defined over end-systole left 
ventricle 
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3. EJECT: Ejection fraction calculation (dynamic method) 



EJECT - ejection fraction calculation (dynamic method) 
The ejection fraction can be calculated from dynamic curves 
stored in a GAMMA-11 save area. The tv?o curves necessary are 
curves of the left ventricle in end-diastole and end-systole. 
The technique for defining the regions of interest is the method 
of Kurtz, Ahnberg, Freed, LaFarge and Treves (to be published) : 

1) the method of creating end-diastole and end-systole images 
is described on page 3, section IB. (Figure 6) 

2) display the end diastole frames, smooth the data and lower 
the upper threshold. (See step 2 in program EJECTS) 

3) enter irregular region of interest (IR) and outline the 
shape of the left ventricle. (Figure 3) 

4) display the end-systole frames, smooth the data and lower 
the upper threshold (see step 5 in program EJECTS) 

5) re-enter IR and outline the end-systole left ventricle. 
(Figure 4) 

6) plot the two curves and then save the curves. (Figure 5) 

7) call program EJECT. 

The program will request the save area number and the curves 
corresponding to end-diastole and end=systole. The end- 
diastole curve is displayed and the user selects pairs of points 
corresponding to end-diastole and end-systole points (Figure 6) . 

fr"h^ f:i"ifar«+- "i rrn ■FT"3r«-t-T one; -Fov cii^r«Ti m,^ i y o-F r^o t -n-f-c: -i c i^t-t n+-ci(^ ^-nr? 

the average ejection fraction for all the pairs is output on 
the terminal. 
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FIGURE 3. ROI defined over end-diastole left 
ventricle 



FIGURE 4. ROI defined over end-systole 
left ventricle 





FIGURE 5. Dynamic curves of the two ROIs 



FIGURE 6. Pick pairs of points (1 diastole, 
1 systole) on ROI curve of L.V. 
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QPQS: Left to right cardiac shunt calculation 



Q PQS - left to right cardiac shunt calculation 
The quantitation of left-to-right shunts can be calculated 
from GAMMA-11 curves stored in a save area according to the 
method of Maltz and Treves. The region of interest (ROI) is 
selected over either lung and free of contamination from the 
heart. The curve is plotted and saved and then the program 
can be executed. The curve of the lung is fit to the gamma 
variate function creating a curve simulating pulmonary flow. 
This curve is subtracted from the original curve to create a 
third curve representing pulmonary recirculation. This curve 
is fit with the gamma variate to simulate premature left- to- 
right recirculation. 

The original curve of the lung, the regression curve, the 
recirculation curve and its regression curve are all stored in 
a save area to allow visual observation of the fit. (See 
Figures 9,10,11) 

The program will request the save area number and the curve to 
be examined and also the number of the save register to store 
the resultant curves. The original lung curve can be 
smoothed and have the background subtracted prior to performing 
the gamma variate fit. In many cases, these two operations will 
improve the statistics of the fit and are to be recommended. 
The point representing 70% of the maxim.um on the downs lope of 

the data curve can be calculated either manually or automatically 
by the program. The peak of the recirculation curve can be 



chosen in a similar manner. 



The program will print the values of the four curves on the 
terminal and also the QP/QS and C2/C1 ratios. 

1) display an image of the lungs with the heart subtracted out 

2) smooth the image (S9) and enter irregular regions of interest 
(IR) . (Figure 7) 

3) draw regions around either one or both of the lungs, and 
plot the curves . (Figure 8) 

4) save the dynamic curves (SDn) . 

5) call program QPOS 

NOTE: Program available in BASIC and FOCAL. 



FIGURE 7. ROIs over left and right lungs 




FIGURE 8. Dynamic curves of the two ROIs 





FIGURE 9. Top- Lung curve range that was fit 
Bottom- Gamma variate fit curve 



FIGUlilE 10 « Top- Recirculation curve 

Bottom- Fit of recirculation curve 




FIGURE 11. Original lung curve overlaid with 
the 2 gamma variate curves . 



READY 

RUNNH 

PROGRAM TO COMPUTE QP/QS 

SOURCE J WHICH SAME REGISTER? 

?J. 

WHICH CURVE? 

?1 

DESTINATION: WHICH SAME REGISTER? 

BACKGROUND ERASURE? 

1»YES 2* mi 
?2 
CURVE SMOOTHING? 

.1. ♦ YES 2 ♦NO 

WHICH LUNG IS REPRESENTED? 

!♦ RIGHT LUNG 2* LEFT LUNG 
?:l. 
MANUAL PICK FOR 70% OF" CURVE ;L? 

It YES 2t NO 

PICK PEAK OF RECIRCULATION CURVE? 

1* YES 2. NO 
"?'? 

MOORE ^ 70-29-00? RACy D? 5/23/74 

CURVE displayed: PICK CURVE DELIMETERS 

LEFT BOUND 

RIGHT BOUND 

CORRECTED CURVE DISPLAYED 

ADVANCE PAPER TO NEW PAGE? TYPE <CR> TO CONTINUE 



MOORE i' 70--29-00y RfhCv Dv 5/23/74 



CURUE 



F"1TTED 



CUF^V 



FIT TED 



396 

1019 

1911 

2666 

3120 

3258 

3024 

2591 

2134 

.f. 6 O Q 

1420 
1337 

J .-1 f\ f-v 
J. a'., li. O 

1130 
1075 
999 

C> .'■ •..} 

■•:■' o i:^ 
/ / .' 

7 4 7 
732 

QPQS =: 



.1. i./ 


•.:> 


/ 


19 


'!•' 


6 


26 


s 




31 


2 


7 


••y '"t 


•1 
.1. 


.1. 


30 

26 





9 

9 


21 




4 


17 


n 





13 


1 


1 


96 


9 




69 


7 




49 


.1. 




33 


9 




23 

•1 C" 




"7 




.1. %J 

10 


1 




66 






42 







O -.J 






'•'l / 

15 




1 

109 
368 
530 

/ "y ft 
O O 7 

736 
769 
742 
698 
6S1 
689 



u 

A 













32 

160 

343 

523 

659 

738 

762 

741 

688 

616 

536 



C2/Ci 



"r -..J * V..' .' /It 



RIGHT LUNG CURVE 



CURUES SAUED IN REGISTER 



5. CARD ; Cardiac output, stroke volume and pulmonary 

transit time calculation 



CARD - cardiac output, stroke volume + pulmonary transit time 
calculation. 

The cardiac analysis program is calculated from dynamic flow 

curves saved in a GAMMA-11 save register. The two regions of 

interest (ROI) that are required are defined over the right and 

left ventricles. The curves are plotted and saved and then the 

program can be executed. 

The area of a curve representing the first transit of the bolus 
is required. The right ventricle curve is displayed and then 
the log of the curve is plotted. The user chooses the bounds 
of a least squares linear fit which is performed on the log 
data displayed (Figure 15) . The exponential of the curve is 
calculated and is combined with the original right ventricle 
curve to eliminate the washout portion of the data. This resultant 
curve represents the first transit and the cardiac output and 
stroke volumn is now computed. (Figure 16) 

The program^ will request the number of the cardiac save area 
and the curve numbers of the right and left ventricles. The 
curves that are created in the program are stored in a save 
area also (Figures 16,17). The blood volumn, equilibrium counts 
and heart rate are also requested in the program. 



The program will print the patient index and the current data, 
then return the cardiac output, stroke volumn, and pulmonary 
transit time. The program returns directly to GAMMA-11 and 
redisplays the last image. 

1) acquire a fast (around 10 frames/sec) dynamic study 

2) display an image of the heart and mark an irregular region 
of interest (IR) over the right ventricle and an IR over 
the left ventricle. (Figures 12,13) 

3) plot the curves and save them in a save area (SDn) . 

(Figure 14) 

4) call program CARD 



i^ 




FIGURE 12. ROI over the right ventricle 



FIGURE 13. ROI over the left ventricle 




FIGURE 14. 



Dynamic curves of the left and 
rightj^lintricles 
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FIGURE 15. 



Pick the linear bounds of the 
log curve 



FIGURE 16. Top- Original right vent, curve 

Bottom- Curve representing 1st trans 




FIGURE 17. 



Log curve and exponential curve 
used during the fit 



8t:B F-ILE NAME— -CARD 



READY 

RUNNH 

CARDIAC ANALYSIS F'ROGRAM 

F'ROGRAM WILL CALCULATE CARDIAC OUTPUT^ STROKE VOLUME^ 
AND PULMONARY TRANSIT TIME 

CARDIAC SAUE AREA •=■• 

':> f\ 

: 7 

RIGHT MENTRICLE CURVEI 

?1 

LEFT VENTRICLE CURVE t 

72 

NEW CURVES SAVED IN AREA: 

?8 

BLOOD VOLUMN (LITERS)? 

?7»0 

EQUILIBRIUM COUNTS (COUNTS/SEC) : 

?150 

HEART rate; 

?75 

LOG CURVE OF RAW DATA 

PICK THE LIMITS FOR THE LINEAR FIT 



CARDIAC ANALYSIS SUMMARY 
22-JAN-76 



patient: HASKINSvCLEMy 23-54-92? RAC> Bf 8/4/75 

CARDIAC OUTPUT ■-'■^ 5*88595 LITERS/MIN 
STROKE VOLUMN == 78*4794 MILLILITERS 
PULMONARY TRANSIT TIME :=:= 4.3 SECONDS 



6. HEAD ; Cerebral blood flow analysis 



HEAD - cerebral blood flow program 

Dynamic brain flow studies can be analyzed quantitatively with 
the program HEAD. The user should define regions of interest 
on either side of the head, then plot and store the curves 
in a save area. 

The program asks the save area number and then whether or not 
the curves should be smoothed and normalized. The first pair 
of curves to be compared are requested and the ratio of the 
amplitudes, ratio of the slopes, ratio of the plateaus, and the 
difference in the times to peak are computed. The program asks 
if there are more curves to be compared and the user must 
respond with "YES", if further analysis is required. 

1) acquire a dynamic brain flow study 

2) display an image of the brain (ADn) , smooth the image (S9) , 
and enter regions of interest (RI) or (IR) 

3) mark regions on either side of the midline (Figure 18) 

4) plot the curves (maximum of 12) and save the cuirves in a 
save area (SD^) (Figures 19,20) 

5) call program HEAD 

NOTE: Program available in BASIC and FOCAL 





FIGURE 18. ROIs marked over each hemisphere 



FIGURE 19. Dynamic curves of the ROIs 




FIGUTRE 20. Dynamic curves after program HEAD 
(smoothed and normalized) 
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OLD FILE NAME-- HEAD 



READY 
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READY 



7. XEBH ; Percent distribution of regional lung ventilation, 

equilibration and perfusion 



XEBH - percent distribution of regional lung ventilation, 
equilibration and perfusion 

Xenon lung studies can be analyzed by using the program 

XEBH. The user should define regions of interest over the 

lungs. It is assumed that there will be an even number of 

ROI's, with the first half of them over the left lung and 

the second half over the right. Normally three regions (top 

of lung, bottom of lung, whole lung) for each lung is 

sufficient, however XEBH can support up to 6 ROI per lung 

(total of 12 curves) . The curves are calculated and saved 

in a save area. (Figure 22) 

The program requests the save area number and then displays 
the last curve calculated. The user is asked to interactively 
choose the ranges for perfusion, equilibration, background 
and ventilation (Figure 23). The program will print the 
patient index and then the percent distribution for each of the 
curves. The V/Q, V/E, and Q/E ratios are also calculated for 
every curve. 

1) acquire a dynamic lung study 

2) display an image of the lungs, enter irregular regions of 
interest (IR) . 

3) mark up to 6 regions over the left lung 

4) mark up to 6 regions over the right lung (the number of 
regions for each lung should be equal) (Figure 21) . 



5) plot the curves and save them in a save area (SDn) (Figure 22) 

6) call program XEBH 

NOTE: Program available in BASIC and K)CAL. 
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FIGURE 21. 6 ROIs marked over the lungs 



FIGURE 22. Dynamic curves of the lungs 




FIGURE 23. Pick the ranges for perfusion, 

ventilajUl^n, equil., arnlBfeackgrd. 
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8. RENO ; Renogram analysis 



RENO - renogram analysis 

Renal studies can be analyzed within the curve analysis program 
in BASIC/RT-11. The user should define 3 regions of interest. 
The first two regions are over the left and right kidneys 
and the third ROI is an area representing the background. 
Normally, the background ROI is chosen in a small area between 
the left and right kidneys (Figure 24) . The curves are 
calculated and stored in a GAMMA-11 save area. 

Run the general purpose curve analysis program CURVE and 
execute the macro RENO (MR RENO) . The program requests the 
save area number and the locations of the three essential 
curves. The 'pick' function is used to mark the bounds of the 
integrals to choose a point on the upslope, to choose a peak, 
and to choose the half life (Figure 26) . 

The program will subtract a normalized background curve from 
the two renograms, then print the total counts, maximum counts 
per minute, time to peak, peak to half-life, slope to peak, 
and slope from peak to half-life for both curves. Left-to- 
right ratios for these values are also returned. Finally, a 
new save area is created containing the three original curves 
and the adjusted left and right renograms with the background 
subtracted out (Figures 27,28). 



1) acquire a dynamic renal study 

2) display an image of the kidneys and enter irregular region 
of interest (IR) . 

3) mark irregular regions over the two kidneys and a third 
region representing the background (Figure 24) . 

4) plot the curves and save them in a save area (SDjj) (Figure 25) . 

5) call program CURVE. 

6) execute the program CURVE and for the 1^"^ command insert: 
MR RENO <CR^. This will execute the macro RENO which is 
set up for renogram analysis. 
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FIGURE 24. 



ROIs over the left and right 
kidneys, and a background ROI 




Fj-GURE 25. Dynamic curves of the 3 ROIs 





FIGURE 26. 



Pick the bounds of the integral, 
then mark the peak, a pt. on the 
upslope, and the half-life. 



FIGURE 27. Left renograms : original curve and 
curve with background subtrabted 




FIGURE 28. Right renogram: original curve 
and adjusted curve 
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9. FILTER; Rib removal of myocardial infarct scintigrams 



FILTER - Rib removed of myocardial infract scintigrams (FORTRAN) 
Imaging of myocardial infracts is subject to distortion 
caused by activity in the ribs overlying the heart. The ribs 
appear as a series of bars running horizontally across the 
image. The program will perform a band-reject digital filter 
on the scintigram and an additional lag filter on every 
column of the image. 

The user must store the image to be filtered in a GAMMA-11 
save area. Mien the program begins, the user inputs the RT-11 
filename of the save area (e.g. SVARJ2fl . SYS , SVAR(2r2 . SYS , etc.). 
The program also requests the centre frequency, the half- 
bandwidth, and the lag constant. The centre frequency is 
determined by counting the number of ribs appearing in the 
scintigrams and extrapolating the number over the entire 
matrix. The half-bandwidth is usually about 30-40% of the 
centre frequency and the lag constant used to smooth the data 
is usually around 50. 

The program is currently designed for 128 x 128 matrices. 
Minor modifications can easily be made for 64 x 64 images, 
but matrices with low resolution do not benefit significantly 
from the filtering technique. 



The digital filter program is an adaption of Stokely's program 
and uses the algorithm of Otnes and McNamee. 

1) Acquire a static image of the myocardium 

2) rotate the image so the ribs are horizontal (ROj^) , (Most 
scintigrams will not have to be rotated) and save the image, 

3) call program FILTER. (CTRL/C 

.R FILTER) 

NOTE: Program in FORTRAN only. 



R FILTER 



IMAGE FILTER PROGRAM 
INPUT SAVE AREA NAME 

*SVAR01.SYS 

PLEASE CHOOSE CENTRE FREQUENCY (:L™64) 
6 

PLEASE CHOOSE HALF BANDWIDTH CI. ••128) 



LAG CONSTANT 
50 

STOP - 



10. FUNC: Functional image program 



FUNC - functional image program (FORTRAN) 

The program FUNC will create a GAMMA-11 patient file consisting 
of a single image which in some way is representative of a 
GAMiVLA-ll dynamic flow study. In this program, either maximum 
counts or time of maximum counts are the criterion for 
creation of the functional image. 

The user inputs the RT-11 filename of the dynamic study and 
the name of the new functional image patient file. (NOTE: 
GAMVLA-ll patient files are required to have extensions beginning 
with an 'X' and followed by two digits 1x00, X01, ...^). The 
program requests which criterion the user wishes for this 
image and the user responds with a 1 for time or a 2 for 
maximum intensity. The functional image can now be examined 
within the GAMIVLA-ll system. 

1) acquire a dynamic study 

2) call program FUNC (CTRL/C 

.R FUNC) (Figures 29,30) 

NOTE: Program in FORTRAN only. The name of the functional 

image file must not be the name of an already existing 
RT-11 file on that device. 




FIGURE 29. Functional Image based on time 




FIGURE 30. Functional Image based on counts 



R FUNC 
INPUT F'ATIENT F'lLE AND NEU FILE NAHE 

*RKi:CRESSW»XOO 
*MAX*XOO 

DO YOU WANT TINE QR INTENSITr FOR THE HATRIX? 
TYPE 1 FOR TIMEj'2 FOR MAX* INTENSITY 



STOP 



R FUNC 

input patient file and new file name 
5|crki:ewingm*xoo 

)|cTIHE*XOO 

DO YOU WANT TIME OR INTENSITY FOR THE MATRIX 
TYPE 1 FOR TIMEf2 FOR MAX* INTENSITY 



STOP 



ADDENDUM 



CURVE Dynamic Curve Analysis 



CURVE is a general purpose BASIC program which allows the user to 
display and manipulate dynamic curves. These curves have been 
saved in a save area and are directly accessible via the GAMMA-11 
BASIC subroutine calls (see Chapter 13) . Addition, subtraction, 
division, and multiplication of curves are possible. Normalizing 
a curve with respect to another, smoothing, plotting the exponetial 
or the log, and integrating are other functions immediately executed 
in this program. Certain functions return a specific value (e.g. 

integration; curve fitting) and these values may be stored in one of 
20 parameters available to the user. CURVE allows the user to 
create macro programs from the list of commands, and also allows 
the user to overlay his own BASIC programs to the nucleus supplied 
by CURVE. 

The first command usually executed during analysis will be to 
read in the save area (RS nl) . Then a dynamic curve should be 
displayed on the terminal (DCnl) . Any curve analysis command 
executed before these two commands are issued, will result in 
unpredictable results and probably fatal errors. 

The other alternative at the start of CURVE is to type the response 
"HELP". This will list all of the available options in the program 
as well as give some general instructions in the operation of the 
system. 

Once the save area and the first curve have been chosen, the user 
can then choose from the list of 33 possible instructions. The 
parameters for these instructions (nl) can be either real values 
or the storage parameters available in CURVE. They are referenced 
as PJ2f,Pl, . . .P9 and Rf2f,Rl, . . .R9. 

Curve Analysis Commands 

1. Input/Output Operations 

a) Read Save Area (RS nl) - read save area nl into the buffer. 
This command must be executed before any analysis can be 
performed on the curves . 

b) Write Save Area (WS nl) - Write the contents of the buffer 
into save area nl. This command stores the results of the 
analysis into one of the nine available dynamic curve save 
areas. This command is usually the last one to be executed 
during analysis. 



c) Write Curve (SC nl) - Write the curve currently on the display 
into the buffer as curve nl. Any curve from Gf to 11 may be 
created or changed in the buffer with this command. Curve 
12 is used by CURVE as the scratch curve. 

2. Display Curves 

a) Display Curve (DC nl) - plot curve nl on the display. This 
command should be executed immediately after reading a save 
area. Then proper analysis of the curves can be guaranteed 

b) Display Overlaid Curve (PO nl) - plot curve nl on the 
display but do not erase the previous curve. The scale for 
the Y- axis remains the same for this curve as it was for the 
previous one. Therefore the overlaid curves may have a 
scale either too large or too small in order to accurately 
view it. There is no limit to the number of overlaid curves 
than can be viewed on the display terminal. 

c) Display Normalized Curve (NO nl) - plot curve nl overlaid 
on the terminal after first normalizing. This function will 
adjust the counts of the curve by the ratio of the number of 
cells of the area currently displayed to the number of cells 
of the curve to be overlaid. 

3. Printing Curves or Parameters 

a) Print Curve (PT) -- print on the output device the curve 
currently being displayed. The scratch curve (12' is drawn 
on the output device and the frame number and the count 
rate are also listed. 

b) Print Parameters (PR nl, n2, . . . ,nl(2f) - print the parameters 
listed on the output device. There may be up to 10 parameters 
listed with this command, and there are a total of 20 storage 
parameters available for the user. The first IGf parameters 
are designated 'pgf ', 'Pi ', 'P2 ',..., 'P9 ' , and the next 10 
parameters are ' RSf' , 'Rl ',,.,,' R9 ' . These parameters may 

be used as either input or output variables in the various 
commands within CURVE. 
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4. Arithmetic Operations 

a) Add Curve (AD nl) -add curve nl to the scratch curve (12) 
and redisplay on the display terminal. The only curve 
actually changed is the scratch curve and the contents of 
the original curves are not altered. 

b'l Add Constant (CA nl) - add the constant nl to the scratch 

curve (12) and redisplay. The appearance of the curve will 
not be altered as a result of this command. 

c) Add Normalized Curve (AN nl) -add curve nl to the scratch 
curve (12^ after first normalizing the curve. Normalization 
is a procedure whereby the counts of curve nl are scaled by 
the ratio of the number of cells of the current curve to 



±1S^II.U^\ 



d) Subtract Curve (SB nl) - subtract curve nl from the scratch 
curve (12i and redisplay. The buffer curve after the 
subtraction will not contain any negative numbers. All 
negative numbers are set to zero. 

e) Subtract Constant (CS nl) - subtract the constant nl from 
the scratch curve (12) and redisplay. Negative values are 
set to zero. 

f) Subtract Normalized Curve (SN nl) - subtract curve nl from 
the scratch curve (12) after first normalizing. Normaliza- 
tion scales the counts of curve nl by the ratio of the number 
of cells of the displayed curve to the number of cells of the 
curve to be subtracted. Negative values are set to zero. 

Multiplicative Operations 

a) Multiply Curves (MT nl) - multiply curve nl by the scratch 
curve (121 and redisplay on the terminal. The scratch curve 
is the only curve altered by the command. 

h) Multiply Constant (CM nl) - multiply the constant nl to 
the scratch curve and redisplay the curve. The appearance 
of the curve will not be altered as a result of this command. 

c) Multiply Normalized Curve (MN nl) - multiply curve nl to the 
scratch curve after first normalizing the curve. Normaliza- 
tion scales the counts of curve nl by the ratio of the number 
of cells of the displayed curve to the number of cells of the 
curve to be multiplied. 
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d) Divide Curves (DV nl) - divide the scratch curve (12) by 
curve nl and redisplay on the scope. If either the dividend 

(curve 12) or the divisor (curve nl) is zero at a particular 
point, the resultant point of the curve is also set to zero. 

e) Divide Constant (CD nl) - divide the scratch curve by the 
constant nl and redisplay. Negative constants are not allowed 

f) Divide Normalized Curve (DN nl) - divide the scratch curve 
by curve nl after it has been normalized. Normalization 
scales the counts of curve nl by the ratio of the number of 
cells of the displayed curve to the number of cells of the 
curve that was indicated. 

6. Curve Manipulations 

a) Plot the Log (LG) - plots the log (base e^ of the scratch 
curve (12) on the scope. This log curve is now the scratch 
curve in the buffer file. 

b) Plot the Exponetial (EP) - plots the exponetial of the scratch 
curve on the scope. It is very easy to exceed the bounds 

of allowable counts with this function. It is recommended 
that the Gommand be issued only when the maximum count of 
the scratch curve is very small. (e.g. after the log is 
displayed, or after the curve has been scaled down via the 
constant divide command) . 

c) Smooth the Curve (SM) - smooth the scratch curve (12) and 
redisplay. The curve is smoothed with a weighting of 1-2-1 
for every point and then is saved in the scratch area. 

d) Integrate the Curve (IN nl,n2,i) - integrates the dynamic 
curve displayed over a specific designated range. The first 
two parameters (nl,n2) are input variables which determine 
the left and right bounds of the integral. The integral is 
computed and returned in the third parameter (i) . The 
parameter i must be one of the storage parameters, since 

it is used solely as an output variable. 

7. Linear Least-Squares pi-t (FT nl,n2,S,Y) 

A linear least-squares fit is performed on the scratch curve 
over a given range of points. The first two parameters (nl,n2) 
are input variables which determine the left and right bounds 
of the fit. The fit is then computed and it returns two values, 
the slope (S) and the Y-intercept (Y) . Therefore, the third 
and fourth param.eters m.ust be two of the storage param^eters, 
since they are output variables that returns a specific value. 
The third parameter will contain the slope of the line computed, 
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and the fourth will contain the Y-intercept. 

8. Point Manipulations 

a^ Pick an X value (PX nl^ - pick a particular X value for the 
curve displayed on the scope. The cursor will appear over 
one of the points of the curve and the user can determine 
which point should be marked. Typing an "L" moves one point 
to the left, typing a "R" moves one point to the right. 
Typing a "J" will move the cursor 10 points in the last 
direction indicated. When the user decides upon a point, 
type an "M" and that point is chosen. The parameter nl 
will return the X value (frame number) of the point in the 
storage parameter designated. 

b'l Pick a Y value (PY nl' - pick a particular Y value for the 
curve displayed on the scope. The cursor will appear over 
the points and the user can determine which point to mark. 
The keyboard commands "L", "R", "J", and "M" move the cursor 
in the sam.e manner as the command PX. lA/hen the user marks 
the given point (by typing "M") , the parameter nl stores the 
Y value (count rate"* of that point. 

c) Change a point (CH nl n2) - changes a point to a new value. 

The point to be altered is determined by the first parameter (nl) 
and the second parameter contains the new value to be stored 
at that point (n2) . 

NOTE: The function PX is frequently used in conjunction with the 
linear fit command (FT) and the integration command (IN) . 
By calling PX twice one can store the left and right bounds 
of the following operations into the area of storage. 



1) PX Pi; PX Rl 
IN Pl,Rl,P2 



2) PX Pi; PX Rl 
FT Pl,Rl,P2,R2 



MACRO Commands 



A macro program is a sequence of curve analysis commands that are 
stored together under one RT-11 filename. The program CURVE can 
then execute all of the commands listed in the macro via one command: 
Macro Run (MR name) . CURVE has all of the editing features necessary 
to create, edit, append, print or execute these programs within its 
system. Once a macro program has been created, it is stored as a 
RT-11 file and can be used as often as necessary during any future 
curve analyses. 
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1. Macro Create (MC name) - creates a macro program from the curve 
analysis commands. All macro commands are illegal within a 
macro. There is no limitation upon the number of commands 
within a macro and the number of commands per line is only 
limited by the number of characters allowed per line (usually 
80 or 132) 

2. Macro Edit (ME name) - edits the macro program designated. 
The program asks the user which line is to be edited. If the 
user types 0, the program exits. The lines of the macro up to 
the edited one are printed, and then the user inserts the new 
line. After one line has been changed, the user can change 
more lines if it is necessary. After exiting the program, the 
remaining lines are printed on the output device and control 
returns to CURVE. 

3. Macro Append (MA name) - adds lines to the macro designated. 
The program prints the macro out to the output device and thert 
requests the additional lines of commands. By typing a null 
line, the program exits and returns control to CURVE. 

4. Macro Print (MP name) - print the macro designated on the output 
device and then return control to CURVE. 

5. Macro Run (MR name) - execute the macro designated. Only curve 
analysis commands and the overlay command are valid within a 
macro program. If an error occurs within the macro, an error 
message will appear and the macro will exit. After completion 
of the macro, control returns to the program CURVE. 

Overlay User Programs (OV name ) 

CURVE allows user programs to be overlaid with the BASIC curve 
analysis system. The overlaid program is executed immediately and 
then control can return back to CURVE. The statement numbers must 
begin at 1000 and be incremented by one. The form of the overlaid 
segment should be that of a subroutine instead of a separate program. 
This means that RETURN rather than STOP should be the last executable 
statement. 

In most overlay segments, the user will access some of the storage 
parameters and system variables within the curve analysis system. 
Certain system variables should not be altered in any overlay 
subroutine (e.g. B9$, 19,09$, ...) . The storage parameters are 
readily available in the overlay by directly accessing the array S. 
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P!2f,Pl, . . . ,P9 are stored in S (gf) , S (1) , . . .S (9) , and R!2f, Rl,...,R9 
are stored in 3(10^), S (11) , . . . ,S (19) . jjtj^gj-gfQjre the overlay can 
be easily used to output the results of curve analysis in the user's 
own unique format. 

e.g. IN 1,20,P5; FT 1,20,R4,R57 OV OUTPUT 

lQi00 PRINT "INTEGRAL OF CURVE ="; S (5) 

l^gTl PRINT "SLOPE OF LINEAR FIT ="; S(14) 

Ij2f02 PRINT "Y- INTERCEPT ="; S(15) 
ief^3 RETURN 

CURVE System Variables 

S(19) storage parameters 

X(128) ,Y(128) curve buffers 

A$(28,2) subroutine name array 

D9 (10) , K9 (10) command parameter arrays 

B1(4),E1(4) curve pointer arrays 

B$ command string 

B9$ command string or macro command string 

09$ overlay name in core 

E$,T$ new overlay name 

C$ command 

C9$ command or macro command 

X9 parameter value 

G parameter number 

18 character position of B$ 

19 character position of B9$ 
K$ position of ";" in B$ 

K9 position of ";" in B9$ 

E4 error flag (E4 = 1, then ERROR) 



PLEASE ENTER COMMAND: HELP 

LIST OF COMMANDS 

MACRO COMMANDS? 
HC NAME MACRO CREATE 
ME NAME MACRO EDIT 
MA NAME MACRO ADD 
MR NAME MACRO RUN 
MP NAME MACRO PRINT 

CURME ANALYSIS COMMANDS J 

OM NAME OVERLAY USER PROGRAM 

RS Nl READ SAME AREA Nl 

WS Nl WRITE SAME AREA Nl 

SC Nl SAME CURME Nl 

DC Nl DISPLAY CURVE Nl 

PO Nl DISPLAY CUR ME Nl ( OVERLAY ) 

NO Nl DISPLAY CUR ME Nl (NORMALIZED) 

PT PRINT CURME ON TERMINAL 

PR NlpN2i' * * ?N10 PRINT STORAGE PARAMETERS 

AD Nl ADD CURME Nl 

CA Nl ADD CONSTANT Nl 

AN Nl ADD NORMALIZED CURME Nl 

SB Nl SUBTRACT CURME Nl 

CS Nl SUBTRACT CONSTANT Nl 

SN Nl SUBTRACT NORMALIZED CURME Nl 

MT Nl MULTIPLY BY CURME Nl 

CM Nl MULTIPLY BY CONSTANT Nl 

MN Nl MULTIPLY BY NORMALIZED CURME Nl 

DM Nl DiiMIDE BY CURME Nl 

CD Nl DIM IDE BY CONSTANT Nl 

DN Nl DIM IDE BY NORMALIZED CURME Nl 

LG PI...OT LOG OF THE CURME (BASE E) 

EP PLOT EKPONTIAL OF THE CURME 

FT NlyN2«N3>-N4 LINEAR LEAS T-SQUARES FIT 

N 1 ■■■■■■■■ 1... !;.. l" T B U N D ? N 2 ■■■■^■- R I G !-i T B U N D ? N 3 :=== S L. P E ? N 4 ■■■^■- Y •- 1 N T E R C E P T 

PX Nl PICK X-MALUE 

PY N:i PICK Y-MALUE 

CH Nli'N2 CMANGE POINT N:l TO N2 

SM SMOOT!"i THE CURME 

IN N 1 i. N 2 V N 3 I n T t G R A T E C U R M E 

N 1 ■-■ l. E f T B U N D y N 2 - R I G H T" BOUND? N 3 =- 1 N T" E G R A L 

THfrRF ARE 20 SIURAGE PARAMETERS FOR THE USERS 

i-*0 V P 1 p r'2 ? P'3 V . * y r"9 y RO y R 1 y R2 ? , * . R9 

THIrlSE hRE USED IN THE CALLING SEQUENCE OF SOME ROUTINES? 

AND MAY B E F R I N f i;: D W :i: T H J H E ' F R N 1 ? N 2 ? . * . ? N 1 " COMMAND 

THE SEMI£:;GL(:iN (?) MAY BE USED TO SEPARATE COMMANDS ON 

THE SAME LINE* NO OTHER CHARACTER IS RECOGNIZED 

AS A COMMAND TERMINATOR EXCEPT FOR A CARRIAGE RETURN. 

SPACES A1<E REQUIRED BETWEEN A MACRO COMMAND 
AND THE MACRO FILE NAME. FOR CURME ANALYSIS 
COMMANDS THEY AR!;: OPTIONAL. 

USER PROGRAMS ARE UMERLAID AT LINE 1000 WITH 
AN INCREMENT OF 1* STORAGE PARAMETER MALUES CAN BE 
STORED AND READ IN fHE MATRIX S(19>. P0?P1?w<.pP9 
ARE STORED IN S ( ) y S ( 1 ) ? . * . ? S ( 9 ) y AND RO ? Rl ? * 1 . v R9 
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EXAMPLE #1: INTEGRATION OF A CURVE 



RUNiMH 

TYPE 'HELr-' FOR A LIST OF COMMANDS 

i ■■ i... e: a s e e h t e r c o h h a n d j 

? R S 9 ? n C 1 — — 



— read save area 9, display curve 1 



i- ' l. E a S !:■ E N J E. R C M n A N IJ t 
? i- • A F ' 1. ? !•■" X !"' 2 ? I N I"" 1. ? I"' 2 M R 1. 

P L. E A S E. E N r I::! ja M H A N li * 

?PR PI mP2?R1 



pick left and rip:ht bounds, 
then integrate 

iDrint the values 



PARAMETER 1 t 
PARAHEFER 2 i 

PARAHETER 11 



43 



!' ■ 1. . i::; Pi S E E N T E R C H M A N B t 



STOP AT LINE 10 
READY 



EXAMPLE #2 : CURVE FITTING, THEN CREATING A MACROS 



lYPE M-IELP" FOR A LIST OF" COMHANDS 



PLEASE ENTER COHHANIK 
?RS8 ? DC2 P NO 1 



F' L E A s f: e n t f: r c: m m a n d t 

TAD J. , 



F L E A S E F: N r E R C H H A N D t 
TLG-iPX P9 JPX P8 

PLEASE ENTER COHHANDI 

tFT F'yyPBj-Ry 9RB 



P i... E A S E E N r E R COMMANDS 
?PR P9j'P8yR9!'R8 



PARAHETER 9 J 
PARAMETER 8 t 

PARAHETER :l.9 t 

PARAMETER 18 * 



21 

* 053261 

6*08133 



read save area 8, display curve 
2, plot normalized curve 1 

add curve 1 



plot the log, pick 2 x values 



linear least-squares fit 



pr-int the values 



PLEASE ENTER COMMAND? 
TMC DEMOl - 



create a macro 



MACRO PROGRAM DEMOl 
INPUT MACRO PROGRAM 
TYPE A NULL LINE TO EXIT 
?RS8n;iC2 0NOliiADl 
TLGsipX P9yPXP8 
? F' T P 9 J' F' 8 9 R 9 y R 8 
'? C' Hi' I"' o „ v.'.' o „ ui' o „ C' o 



LEASE ENTER COMMAND 



STOP AT LINE 10 
READY 



EXAMPLE #3 : MACRO COMMANDS 



RUNNH 

TYPE 'HELP' FOFi: A LIST OF COMMANDS 



PLEASE ENTER COMMAND? 
?ME DEMOl — — — 



edit macro demol 



MACRO PROGRAM DEMOii. 

EDIT MACROI (TYPE TO EXIT) 

WHICH LINE ? 

?3 ~ 



LINE 1 RS8?DC2?N01 

LINE 2 LG5PX P9yPX 

LINE 3 FT P9!'P8?R9 

INSERT NEW LINE I 
TIN P9n-'8yP7?FT P9 5'P8?R8? 
WHICH LINE ? 



yADl 
P8 
9 R8 

R9 — 



replace line 3 



new line inse±?ted 



TO 
LINE 4 



PR P9yP8?R9)'R8 



PLEASE ENTER COMMAND 
TMR DEMOl 



execute the macro 



MACRO PROGRAM DEMOl 
PARAMETER 9 I 12 

PARAMETER 8 X 29 

PARAMETER 19 % 6*50975 

PARAMETER 18 % * 02571 11 



PLEASE ENTER COMMAND I 

TPR ?7 

PARAMETER 7 J 126*663 

PLEASE ENTER COMMAND? 

TMP DEMOl 



print param.eter 7 
(integral of the curve) 

print the macro 



MACRO PROGRAM DEMOl 

LINE 1 RS8n:iC2iSN01?ADl 

LINE 2 LG?PX P9?PXP8 

L I N E 3 IN !"' 9 ? F' 8 '>?"/? F T I-' 9 

LINE 4 PR P9?P8?R9yR8 






PLEASE ENTER COMMAND 



STOP AT LINE 10 



READY 



EXAMPLE #4 : PRINTING THE CURVE ON THE TERMINAL 



?Vffi'"'''"HELF-'^ FOR A LIST OF COHHANDS 



PLEASE ENTER COMHAND 
?RS9yDC1. ?PT — 



read save area 9; display curve 
iDrint curve on the terminal 
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EXAMPLE #5 : USE OF USER OVERLAY PROGRAM 



RS9?DC2 



read save area 9, display curve 2 



P I. E A S ii;: E H T E R C; H M A N B i 

i^PX P6^PX P7DIN P6?P7?R2yOU IiEMu2 -J — pick left and right bounds of the 

integral, intesrr ate, call user 
LEFT BOUNDS OF INTEGRAL ■■- 6 overlav -oroerram ! 

RIGHT BOUNDS OF INTEGRAL "" 40 

PATIENT NAMEJCAYERy 49-"06-8;L y CARDIAC? D? 11/15/73 
INTEGRAL ■■"■ 5245 

PLEASE ENTER COMMAND J 

'T' 

SrOP AT LINE 10 

READY 

OLD 

G;. D FILE NAME--DEM02 

KEADY ■ " list the user overlay program 

LJSTNH 

tOOO PRINT \PRINT "LEFT BOUNDS OF INTEGRAL ■■^^'"?S(6) 

1 i !■■' R I N T " R 1 G H T B U N D S F' I N T E G R A L "■•■ " ? S ( 7 ) 

1002 CALL " GSAG " ( 1 ? A3$ ) 

1003 PRINT XPRINT "PATIENT NAMErSA3$ 
1 4 F R I N T ' 1 N T E G R A L -- " y S (12) 

100:5 RETURN 

READY 



